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A static microscopic study of potential-energy surfaces within the Skyrme-Hartree-Fock-plus-BCS 
model is carried out for the ^^^Fm and ^^^Fm isotopes with the goal of deducing some properties 
of spontaneous fission. The calculated fission modes are found to be in agreement with the experi- 
mentaly observed asymmetric-to-symmetric transition in the fragment- mass distributions and with 
the high- and low-total-kinetic-energy modes experimentally observed in ■^^^Fm. Most of the results 
are similar to those obtained in macroscopic-microscopic models as well as in recent Hartree-Fock- 
Bogolyubov calculations with the Gogny interaction, with a few differences in their interpretations. 
In particular an alternative explanation is proposed for the low-energy fission mode of ^^*Fm. 

PACS numbers: 21.60.Jz,24.75.+i,27.90.+b 



I. INTRODUCTION 

Before the work of Brandt and collaborators ^ in 
1963, very little was experimentally known about de- 
tailed spontaneous-fission properties of isotopes other 
than ^^^Cf. These authors measured the mass- yield and 
kinetic-energy distributions for the spontaneous fission 
of ^^^Fm and obtained a well-marked asymmetric mass- 
yield curve. Later on, Balagna et al. Ql and John et al. p| 
showed independently that the experimental fragment- 
mass distribution for the spontaneous fission of ^^^Fm is 
essentially symmetric, with a very slight dip in a broad 
peak centered at the fragment mass 127-128. John and 
collaborators 3] also investigated the thermal-neutron- 
induced fission of ^^^Fm and obtained a sharper symmet- 
ric pattern. Based on the results from the Argonne group 
published the following year [4l| and showing the asym- 
metric character of the spontaneous fission of ■^^^Fm, they 
concluded that "symmetric mass division in low-energy 
fission of heavy actinides appears abruptly at ^^^Fm". 
This was confirmed later in the measurements by Hulet 
and collaborators [ESQ ^^'^ found very narrow sym- 
metric mass distributions in the spontaneous fission of 
258Fm and ^s^Fm. 

Another remarkable spontaneous-fission property that 
rapidly changes along the Fm isotopic chain is the frag- 
ment kinetic-energy distribution: it is very well repro- 
duced by a single Gaussian for A ^ 256 with an av- 
erage total kinetic energy TKE that follows the Viola 
systematics [y, |8j , whereas heavier isotopes have a much 
higher TKE- value. This is particularly the case for ^^^Fm 
where Hulet et al. j3] found that the total-kinetic-energy 
distribution has a non-Gaussian shape which can be un- 
folded into two Gaussians. These authors named this 
behavior "bimodal fission" and showed that it corre- 
sponds to two different fission modes: one is a high-energy 
mode (TKE = 230 MeV) associated with a narrow sym- 
metric mass distribution, the other one is a low-energy 
(TKE = 205 MeV) form of fission with a much broader 
(still symmetric) mass distribution, which even reverts to 



asymmetric when spontaneous-fission events associated 
with lower energies (TKE ^ 200 MeV) are selected. In- 
terestingly both modes have about the same abundance. 

From the theory side, a number of authors investi- 
gated these two spontaneous-fission properties in Fm iso- 
topes. On the one hand, Lustig, Maruhn and Greiner \\Sj\ 
studied the transition in the mass-distribution patterns. 
They calculated the mass-yield curves of even Fm iso- 
topes from A = 254 to A = 260 and reproduced the 
observed transition. On the other hand, several groups 
investigated the bimodal fission of ^^^Fm, from the char- 
acteristics of its potential-energy surface. Four of them 
made use of a macroscopic-microscopic model, relying on 
a liquid-drop contribution to the binding energy (with dif- 
ferent nuclear surface parametrizations) and the Strutin- 
sky method with various single-particle potentials to cal- 
culate the microscopic correction to the macroscopic en- 
ergy Hl 12, JJi 14; I5J- Even though all of them did 
not agree on the height of the outer saddle point relative 
to the ground state of the fissioning nucleus, they all ob- 
tained two fission valleys leading to two different families 
of prescission shapes, namely elongated slightly asymmet- 
rical configurations, and compact symmetrical configura- 
tions corresponding to two nearly spherical '^^^Sn nascent 
fragments. They were dubbed respectively as the "old 
path" and the "new path" in Refs. J12. Il3 |. The exis- 
tence of these two fission valleys is in agreement with the 
experimental observation of two energy modes. Pashke- 
vich [iJI even showed that the difference between the 
Coulomb interaction energy between the nascent frag- 
ments estimated in both valleys for different neck radii 
(about 30 MeV) was reasonably close to the TKE differ- 
ence between the fission modes (25 MeV). 

More recently, two groups of authors performed cal- 
culations of potential-energy surfaces, covering the re- 
gion of saddle points and beyond, within microscopic self- 
consistent mean-field approaches including residual pair- 
ing correlations and using density-dependent phenomeno- 
logical effective nucleon-nucleon interactions. On the one 
hand, within the Hartree-Fock-Bogolyubov (HFB) ap- 
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proximation with the DIS Gogny interaction, Warda and 
collaborators [1^, llS uM carried out an extensive study 
addressing both spontaneous-fission properties discussed 
in the first two paragraphs. However some of their argu- 
ments do not seem to be consistent with their results, as 
will be shown in Sect. 3. On the other hand, Staszczak 
and coworkers 19] implemented the Hartree-Fock-plus- 
BCS (HFBCS) approximation, with the SLy4 Skyrme 
interaction in the mean-field channel and the seniority 
(constant-G) force in the pairing channel. They also ob- 
tained a reflection-asymmetric path and a symmetric one 
in their fission-barrier calculations for the even Fm iso- 
topes from A = 242 to A = 264, but their claimed over- 
all qualitative agreement with existing experimental data 
does not seem consistent with the results presented. 

Finally it is worth mentioning that Asano and cowork- 
ers [23 recently performed dynamical calculations of frag- 
ment kinetic-energy and mass distributions, based on the 
multi-dimensional Langevin equation and a macroscopic- 
microscopic potential-energy surface, for the ^^^Fm iso- 
tope at a compound nucleus excitation energy of 10 MeV. 
They also calculated the mass distributions of ^^^Fm and 
258pj^ at the same excitation energy and obtained a good 
agreement with the experimental data for the thermal- 
neutron-induced fission of ^^^Fm |2J and ^^"^Fm Q. 

Within the HFBCS microscopic approach, the present 
study aims at clarifying several points related to the 
above spontaneous-fission properties of ^^^Fm and ^^^Fm 
from a careful, static study of their potential-energy sur- 
faces, extending the brief discussion in Ref. |22|. After a 
short description of the HFBCS model and the method 
used to explore the energy surface, I will present the re- 
sults in Sect. 3 and discuss them in Sect. 4. Finally I will 
draw the main conclusions in Sect. 5. 



II. MODEL AND METHOD FOR EXPLORING 
THE POTENTIAL-ENERGY SURFACE 



The microscopic approach followed in this study has al- 
ready been applied to calculations of actinide fission bar- 
riers and presented in detail in Ref. [23 ■ It is based on 
the Hartree-Fock-plus-BCS approximation implemented 
with a density-dependent nucleon-nucleon effective inter- 
action. Such an approach is self-consistent owing to both 
the Hartree-Fock approximation itself and the density- 
dependence of the effective interaction. However, since 
very large elongations are involved up to and beyond scis- 
sion (defined for example as in the recent works |2J| and 
[23), the approximate correction term for the zero-point 
rotational motion added in Ref. ^3 does not seem to be 
well suited for such deformations and is not taken into 
account. As for the spurious contribution from center- 
of-mass vibrations to the energy, thoroughly analyzed by 
Bender et al. in Ref. I2g|, only the traditional one-body 
contribution E^ (leading simply to a renormalization 



factor 1 — 1/A in the kinetic energy) is considered since 
the SkM* force parameters were fitted within this frame- 
work. As we can see in Tab.U adding the two-body con- 
tribution E^ ' (perturbatively calculated) to the Hamil- 
tonian changes the deformation energy with respect to its 
ground-state value by about 1 MeV. As long as we deal 



TABLE I: One-body E^^' and two-body Ej^' contributions to 
the center-of-mass kinetic energy in MeV per nucleon at the 
ground state (Q20 ~ 32.5 b) and far beyond the outer saddle 
point in the fission valley (Q20 ~ 400 b) of ^^^Fm (see Sect. 3 
for details). 



Q20 (b) 


eH^/a 


E'^^/A 


total 


32.5 

400.7 


18.72 
18.34 


-13.21 
-12.16 


5.51 
6.18 



with one-cluster nuclear shapes, this effect can be con- 
sidered to be small. However, the center-of-mass correc- 
tion is expected to be important at and beyond scission, 
and a tractable correction appropriate for any configu- 
ration ranging from slightly deformed to well-separated 
fragment shapes has not yet been proposed. 

To illustrate this, let us consider a configuration with 
two identical, well-separated fragments. On the one 
hand, we would expect the energy of the total system 
to be 



E = Ei+E2 + -B^'"'^ 



(1) 



,(int) 



where E^ is the Coulomb interaction energy between 
the two fragments (which expression, not needed here, 
will be given after this discussion). In Eq. ^, Ei denotes 
the energy of fragment i obtained by subtracting the one- 
and two-body parts of the center-of-mass kinetic energy 
of fragment i from the expectation value {H){A/2) of the 
corresponding Hamiltonian: 



E, 



W(A/2)-(^^ + ^i^j. (2) 



Using the following estimates for E^ ' j A and Et^ ' jA as- 
sumed to be independent of A and the deformation co- 
ordinate (expressed as the quadrupole moment Qiq) and 
deduced from Tab. Q] 



£;^^V-4 « 18.5 MeV 
-0.7B^^V^ 



,-.(2) 



E\,^^IA 



we arrive at 



^«2(iJ)(A/2)+B, 



(int) 



11 MeV. 



(3) 
(4) 

(5) 



On the other hand, a common application of the center- 
of-mass correction (even including the two-body part) 



leads to 



E' = {H){A) 



p(i) 



E\^\A) , E)l\A) 



(2)/ 



A 



A 



r(int). 



with {H){A) w 2 {H){A/2) + E^^"''; hence 

^' w 2 {H){A/2) + ^^'"*' - 5.5 MeV. 



(6) 



(7) 



The difference between E, the intuitively correct result, 
and E' amounts to about 5 to 6 MeV in absolute value. In 
this work, where the two-body part of the center-of-mass 
correction is not taken into account, the energy actually 
calculated is 



E" = {H){A) 



eII\a) 



A 
2{H){A/2)+E, 



(int) 



18.5 MeV, 



(8) 



which underestimates E by about 7 to 8 MeV. That 
means that, in the absence of a correct treatment of 
the spurious center-of-mass motion, the approximate cor- 
rected energy calculated here is off by at least 5 MeV at 
and beyond scission, even if the two-body contribution 
is taken into account. Nevertheless for the purpose of 
the present study this is not a serious limitation since it 
does not affect the existence of stable solutions at a given 
elongation. 

Let us now be precise about how the Coulomb energy 
is calculated in this work and what expression of E^ is 
used to estimate the fission-fragment total kinetic energy 



in Sect. 4. Firstly the direct Coulomb energy E^^ 
the fissioning nucleus 



(dir) 



of 



E. 



(dir) 



73„ / j3 /Pp(i")Pp(r') 



d^r / d^r' 



r - r' 



y f dh f d^r'\r-r'\Appiv') 



(9) 



(10) 



where e^ « 1.44 MeVfm, Pp{r) is the (local) proton den- 
sity and A is the Laplacian operator, is computed ex- 
actly by numerical integration of the well-behaved ex- 
pression (|10|) as proposed by Vautherin |23 . Then the ex- 
change part Eq^ ' is calculated by combining the Slater 
approximation |28j| with a kind of local density approxi- 
mation, which leads to the expression 29, 30] 



e: 



(exch) 



e 
3eV3 

4 \ TT 



d^r I d^r'^i^^ 



1/3 



-^ Z dh[p,{r)] 



,4/3 



(11) 

(12) 



where pp(r, r') is the nonlocal proton density (off- 
diagonal matrix elements of the one-body density opera- 
tor). The two densities Pp(r) and pp{r,r'), including the 



BCS occupation probabilities v^ as in Ref. ^27], take the 
form 



Pp 



w = 2E-'El^rnr) 



H/^^|2 



(13) 



i>0 



p,(r,r') = 2 E".'E(^l'^^W)Vi''^(r') ' (14) 



i>0 cr.cr' 



where the sums run over all the pairs of time-reversed 
conjugate states and (pi^ (r) stands for the component of 
spin a of the single-particle wave function associated with 
the single-particle state \i). The Coulomb interaction en- 
ergy between two complementary pieces Vi and V2 of the 
fissioning nucleus is defined by 



4'"*) - e2 



Vl JV2 



Pp(ri)pp(r2) - |pp(ri,r2)| 



ri-ra 



(15) 
and can be decomposed into direct and exchange parts 

4 



n(int dir) n 7-^(int exch) .. ^ 

' ^ ' and E\ ' , respectively: 



E, 



(int) 



E, 



(int dii 



E, 



(int cxch) 
C 



with 



E, 



(int dir) 
C 



/ d^ri pp(ri) / d^v-2, 



e: 



(int cxch) 

c — 



Vl 



Pp(r2) 
ri -r2| 
2 



dhi / d^r2 



V2 



|/3p(ri,''2) 

|ri -r2| 



(16) 

(17) 
(18) 



Eqs. (|14|l and p8l) show that the exchange contribution 
vanishes when the single-particle wave functions are lo- 
calized in either subset Vi or V2, which is the case in 
separated-fragment configurations (at and beyond scis- 
sion), but not before the neck ruptures. Since the ex- 
change part of the Coulomb energy is calculated with the 
approximate expression H12|l . it gives no contribution to 



the Coulomb interaction energy. Indeed E^ 
written as 



(cxch) 



E, 



(cxch) 



E. 



(cxch 1) 
C 



E. 



(cxch 2) 
C 



where 4 is given by 



E. 



(cxch i) 



3e^ 



1/3 



dh [pp{r] 



14/3 



can be 



(19) 



(20) 



Vi 



and is interpreted as the exchange term of the Coulomb 
self-energy of the subset Vi. In consequence, 4™ is 
overestimated when its exchange contribution is calcu- 
lated with the approximation H12() . but by less and less 
as we approach scission. 

In practice the phenomenological Skyrme SkM* inter- 
action in the mean-field channel and the seniority force 
in the pairing channel are chosen, with the same set 



of parameters as in Ref. |23|. The Hartree-Fock equa- 
tions are solved by expansion of the single-particle states 
in a cylindrical harmonic-oscillator basis, which needs 
to be appropriately truncated for practical calculations. 
Appropriately means here that, given a truncation pre- 
scription and a maximal size parameter A^o (see, e.g, 
Ref. |3ll|l. the basis parameters for a given deformation 
should be chosen so as to minimize the energy. Following 
Ref. [Sj, I introduce a deformation parameter q — uj^juj^ 
and the spherical-equivalent harmonic-oscillator constant 
h — y^muj/h, where w is related to the oscillator frequen- 
cies (jJz (in the z direction) and uj± (in the plane perpen- 
dicular to the z-axis) through tu^ — lu^ uj'j_ . In particular 
17 oscillator shells {Nq = 16) are included throughout 
this work unless otherwise specified, as this basis size 
was shown in Ref. J23 to be sufficient to within about 
0.1 MeV for the relative energies (with respect to the 
ground state) up to the outer saddle point of ^^^Cf. 

Given the variational character of the HFBCS approx- 
imation and the difficulty of computing energy on a mesh 
in an A'^-dimensional deformation space with A'^ > 2, I re- 
sort to implementing the method described in Ref. [221 
where it was applied to calculations of fission paths for 
^"Se. Assuming axial symmetry, a limited number of 
shape degrees of freedom are retained, namely the elon- 
gation of the fissioning system expressed as the axial 
quadrupole moment Q20 or the center-of-mass distance 
-Dc.m. between the (pre-)fragments, the mass asymme- 
try through either the axial octupole moment Qao or the 
heavy (pre-) fragment mass Ah^ and the neck coordinate 
Qn (introduced by Berger et at ^3 and used by Warda 
et al. [l6J)- Under this assumption, the center of mass 
of the fissioning nucleus, located on the symmetry axis 
chosen to be the z-axis, is fixed at the origin of the ref- 
erence frame by adding to the Hamiltonian a constraint 
on the expectation value of z. The definitions of the re- 
tained shape coordinates can be found in Ref. _32l| , except 
for Ah 



Ah = max{ Aright, ^loft} 



(21) 



where the mass of the right and left fragments are respec- 
tively defined by 



Aight = / rf r p(r) 

^Icft —A — Aright • 



(22) 
(23) 



In Eqs. (j^ and if^ . p(r) denotes the nuclear density 
(neutron plus proton contributions) and A the mass of the 
fissioning nucleus. The neck abscissa Znock is defined here 
as the z-value at which the nuclear density integrated 
in the perpendicular plane is minimal. This definition 
holds only for sufficiently necked-in nuclear shapes and 
corresponds to the value of z for which the neck radius 
is minimum. The light-fragment mass A^ is thus simply 
given by Al = A — Ah- The number of protons in the 



heavy and light fragments Zh and Zl are defined in the 
same way. 

Interpreting the fission modes in terms of valleys of the 
potential-energy surface, we do not need to explore the 
whole surface to find these valleys and a constrained vari- 
ational approach like the HFBCS approximation seems 
suitable to this purpose. However, two limitations should 
be kept in mind. On the one hand, a partial exploration 
of the energy surface might not lead to the lowest valleys 
at a given elongation (the driving coordinate here), on the 
other hand such an approach does not guarantee to find 
the lowest saddle point between a pair of local minima 
(in a given deformation space), as already pointed out 
in Ref. 32]. It will therefore be verified a posteriori, by 
comparison with experiment, that the valleys obtained 
are the most relevant ones. An important point when 
using this approach is to take great care to check the 
stability of the solutions when looking for the relevant 
valleys. More precisely these solutions should correspond 
to significantly deep local minima for each constrained 
Q20-value in the directions of mass asymmetry and neck 
coordinate. Finally, the point in a given fission valley 
where the solution corresponds to the larger constrained 
elongation is called the exit point of the fission valley. 
One can similarly define the exit point of a fusion val- 
ley upon considering the solution that corresponds to the 
smaller constrained elongation. 



III. RESULTS 

The model and method presented in the previous sec- 
tion are applied to the ^^^Fm and ^^^Fm isotopes. Earlier 
calculations [23i | were performed for these heavy nuclei in 
the same framework except that left-right symmetry was 
imposed (in addition to axial symmetry). The results 
reported here are obtained by releasing this constraint. 

For the ^^^Fm isotope, only left-right reflection sym- 
metric solutions are found between the spherical shape 
at Q20 — and Q20 ~ 160 b. The corresponding 
deformation-energy curve is plotted in Fig. Q as a solid 
line. A superdeformed minimum lies just beyond the in- 
ner fission barrier, around which preliminary calculations 
seem to indicate some softness in the direction of left-right 
asymmetric deformations (with a possible local minimum 
at a finite Qao-value). Beyond this minimum, two kinds 
of valleys exist: the so-called fission valley, correspond- 
ing to one-body-shaped configurations, and the so-called 
fusion valley, corresponding to two separated fragments 
(see, for example, Refs. |33ll34|). 

On the one hand, only one fission valley is obtained 
(dashed line labeled "asym. EF" - asymmetric elon- 
gated fission - in Fig. ^|, along which the ^^^Fm nu- 
cleus exhibits left-right asymmetric shapes. The valley 
stretches from Q20 ~ 140 b to Q20 ~ 420 b, where the 
one-body-shaped solution becomes unstable against neck 
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FIG. 1: Potential energy of deformation _Edcf as a function of 
the quadrupole moment Q20 along the different fission paths 
obtained for the ^ Fm isotope. The solid line from to 160 
barns (b) represent left-right symmetric solutions including 
the ground-state (GS) and the superdeformed (SD) minima. 
The solid and the dotted lines beyond 250 b correspond to 
symmetric and asymmetric compact fission (CF) paths, re- 
spectively. The dashed line is the bottom of the asymmetric 
elongated fission (EF) valley. 



rupture. As we can see in Fig. |21 showing slices of the 
potential-energy surface in the Q^q direction at three dif- 
ferent fixed elongations Q20, the transition between the 
symmetric path from the superdeformed minimum and 
the asymmetric fission valley occurs smoothly around 
Q20 = 160 b. At the exit point of the fission valley, 
the nascent fragments are calculated to be approximately 
mXe and ^i^Pd. 

On the other hand, two fusion valleys corresponding 
to the symmetric ^go^n-f-^ggSn and sHghtly asymmetric 
"'^4gIn-|-^|5Sb fragmentations are found in the potential- 
energy surface. They are represented in Fig. Q] as a solid 
line and a dotted line, labeled "sym. CF" (symmetric 
compact fission) and "asym. CF" (asymmetric compact 
fission), respectively. Although these valleys are plotted 
only for elongations less than Q20 = 375 b, they exist for 
larger elongations. In contrast their upper ends represent 
exit points as defined in the end of the previous section. 
The two fusion valleys can also be visualized in Fig. 
showing a cut in the heavy-fragment mass Ah direction 
at a fixed elongation Q20 = 280 b chosen as an example. 
In this figure, the solid line represents solutions obtained 
by increasing the constraint on Ah from 128 (symmet- 
ric solution) to about 133, beyond which these solutions 
become unstable. Similarly the dashed line corresponds 
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FIG. 2: Cuts in the potential-energy surface along the axial 
octupole moment Qzo, direction at different fixed elongations 
(320 for ^^®Fm. The full circles represent calculated points, 
whereas the solid lines are drawn as eye guides. 



to solutions obtained by decreasing the constraint on Ah 
from about 141 to about 131, below which this kind of 
solutions does not exist. The two minima at Ah = 128 
and Ah = 133.5 are therefore not connected through a 
continuous curve in Fig.|21 from which we learn that one 
or several additional degrees of freedom are missing in the 
description of this region of the energy surface (for exam- 
ple the fragment deformations). Nevertheless this does 
not cast any doubt on the existence of the two fusion 
valleys. Regarding the exact location of the asymmet- 
ric minimum on the dashed curve of Fig. 13 the present 
HFBCS calculations indicate that the mass of the heavy 
fragment varies very little in the asymmetric fusion valley, 
with an average integer value of 133. 

Let us now turn to the ^^®Fm isotope. Its potential- 
energy surface presents some features in common with the 
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FIG. 3: Cut of the potential-energy surface along the Ah 
direction at a fixed elongation Q20 ~ 280 b for ^^^Fm. The 
solid line represents solutions obtained by increasing Ah from 
128 (symmetric solution) and the dashed line is obtained by 
decreasing Ah from about 141. 



FIG. 4: Same as Fig.Qfor the ^''^^Fm isotope. FuU circles and 
triangles are calculated points along the symmetric compact 
fission ("sym. CF") and symmetric elongated fission ("sym. 
EF") paths, respectively. 



one of ^^^Fm from the spherical point to the also present 
superdeformed minimum, with only symmetric solutions 
- corresponding to the solid line in Fig. 21 ~ in this elon- 
gation range. As in ^^^Fm, octupole softness is observed 
around the top of the inner fission barrier. An asymmet- 
ric fission valley emerging in the vicinity of the superde- 
formed minimum is also present in ^^®Fm and plotted as a 
dashed line in Fig. 0] The nascent fragmentation around 
the exit point, namely approximately "'^|4Xe+^4gPd, is 
similar to that obtained in ^^^Fm. However there are 
two major differences between ^^^Fm and ^^^Fm. The 
symmetrical solutions in ^^^Fm not only remain stable 
against left-right asymmetric deformations until scission, 
but they also give rise to two different families of nuclear 
shapes. The solutions with compact shapes constitute 
the symmetric compact fission (CF) path marked with 
full circles in Fig. 21 (labeled "sym. CF") and those asso- 
ciated with elongated configurations form the symmetric 
elongated fission (EF) path marked with full triangles 
(labeled "sym. EF"). 

A word of caution should be said here about the sym- 
metric EF path. As mentioned in Sect. 2, the present 
calculations are done using a harmonic-oscillator basis 
size A'o — 16 for the expansion of single-particle states. 
Even though this is large enough to show the existence of 
the symmetric EF path (that is, its stability with respect 
to left-right asymmetric deformations), it is not the case 
when addressing quantitative questions like the location 



of the exit point. This is why axial and left-right sym- 
metric calculations are performed with a larger basis size 
A'o — 20. It then becomes possible to scan the potential- 
energy surface on a mesh (I?c.m., Q2f), where Dc.m. and 
Q2{ stand, respectively, for the center-of-mass distance 
and the fragment axial quadrupole moment (which is of 
course the same for each fragment because of the reflec- 
tion symmetry) . Since we deal here with very large elon- 
gations of the fissioning nucleus, the centers of mass of 
the pre-fragments (and a fortiori of the separated frag- 
ments) are well defined and it becomes more physical and 
intuitive to consider -Dc.m. rather than Q2o- It is worth 
adding here that Q20^ Q21 and Dc.m. are not independent 
of each other since they obey the relation 



Q20^2Q2i + -Dl, 



(24) 



where A is the mass number of the fissioning nucleus. 

In practice, it is extremely time consuming to opti- 
mize the basis parameters b and q at each grid point. 
Therefore I resort to an approximate procedure in which 
b is optimized using a smaller basis corresponding to 
Nq = 16 at an elongation Q20 ~ 190 b somewhat smaller 
than the elongation at which the symmetric EF valley 
appears. The resulting value of b is then used over 
the whole range of deformations covered by the mesh 
(10 fm ^ Dc.m. sC 23 fm and -30 b s$ Qaf < 60 b). As 
we can see in Tab. (^ 0.42 is the approximate optimal 
value for b. As for q, to which the results are not very 



TABLE II: Values of the deformation energy on a {b, q) mesh 
{b in fm"^) at Q20 = 190 b for ^^*Fm, calculated with 17 
major shells (A^o = 16). 
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1.6 


1.8 


2.0 


2.2 


0.38 


-1903.03 


-1903.28 


-1903.75 


-1903.67 


0.40 


-1903.85 


-1904.04 


-1904.21 


-1903.93 


0.42 


-1903.99 


-1904.05 


-1904.21 


-1903.66 


0.44 


-1903.61 


-1903.66 


-1903.91 


-1903.37 


0.46 


-1903.15 


-1903.31 


-1903.56 


-1903.16 



Fm N(,=20 b=0.42 




D. ^ (f m) 



Q2f (b) 



'c.m.v 



sensitive (see Tab. lnj) . especially with the enlarged basis, 
an approximate variation with Uc.m. is taken into account 
(the actual values of q appear in Fig. [3 discussed below). 
The resulting numerical uncertainties are not expected 
to drastically affect the relative position of the symmet- 
ric EF and CF valleys or the nuclear shape at the exit 
point of the EF valley. 

The potential-energy surface obtained for ^^*Fm is dis- 
played in three different forms: as a three-dimensional 
surface in Fig. [S] as a two-dimensional contour diagram 
in Fig. and as a series of cuts at fixed Z?c.m.- values in 
Fig. [3 The deformation-energy curves as functions of 
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FIG. 5: Surface of deformation energy as a function of the 
center-of-mass distance Dc.m. and the fragment elongation Q2t 
for 25*Fm. 



FIG. 6: Contour diagram of the deformation energy as a func- 
tion of the center-of-mass distance Dc.m. and the fragment 
elongation Q2{ for ■^^*Fm. The energy interval between the 
contour lines (dEdcf) is 3 MeV. 



the symmetric EF valley. The deformation-energy curves 
are all continuous except for I?c.in. ^ 21 fm (for which 
only the calculated points are plotted) . Indeed two kinds 
of configurations coexist for 30 b ^ Q21 ^ 40 b. They 
are characterized by different Qw -values: Qn ~ for 
the steep increasing branch (separated fragments) and 
Qn ~ 6 for the other one (very elongated one-body- 
shaped configuration). From Fig.Qwe can approximately 
localize the exit point of the symmetric EF valley at about 
I?c.m. = 22 fm, with Qaf ~ 45 b. 

From these results, we can deduce the variation of the 
deformation energy along the bottom of each valley as a 
function of the driving coordinate Dc.m. (see the upper 
panel of Fig. (S)). It is also instructive to plot the same 
curves as functions of the quadrupole moment Q20 (see 
the lower panel of Fig. IHJ. The variation of the energy is 
of course not affected by the choice of the driving coor- 
dinate since -Dc.m. is a monotonically increasing function 
of Q20 as shown in Fig. El but the relative position of the 
two valleys can differ because of a projection effect. In- 
deed projecting the multidimensional energy surface onto 
a deformation subspace distorts it and the resulting pat- 
tern generally depends on the actual subspace. However, 
we can neglect this effect when comparing two valleys if 
the distortion is weak or the energy difference between 
the two valleys is large. 



Q21 in Fig. Ul exhibit two minima for Dc.m. ^ 12.5 fm. 
One at Q2f ~ is associated with a configuration having 
very few nucleons in the neck (Qn ~ 0), corresponding 
thus to two separated, identical and nearly spherical frag- 
ments (symmetric CF valley) . The other minimum varies 
with Dc.m. from 18 b to 45 b and the associated configura- 
tion has a finite neck radius {Qn ^6): it corresponds to 



IV. DISCUSSION 

A. Asymmetric-to-symmetric transition in the 
fragment-mass distribution 

From the features of the potential-energy surface we 
can obtain some information about the fragment-mass 
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FIG. 7; Cuts in the potential-energy surface along the Q2i direction at different fixed center-of-mass distances Dc.m. for ^^**Fm. 
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FIG. 9: Variation of Dc.m. along both symmetric fission val- 
leys as a function of Q20, obtained with A'^o = 20 and b = 0.42. 



FIG. 8: Deformation energy along both symmetric fission val- 
leys as a function of Dc.m. (upper panel) and as a function of 
Q20 (lower panel), obtained with A^o = 20 and b — 0.42. 



distribution in the spontaneous fission of ^^^Fm and 

258Fm. 

Starting from the ground state of ^^^Fm, the lowest 
and only continuous path leading to scission is asymmet- 
ric in its late stages, where the nascent fragments have a 
fairly constant mass ratio of Ah/Al ~ 140/116 beyond 
Q20 ~ 350 b. Under the assumption that the most prob- 
able fragmentation experimentally observed corresponds 
to the configuration just before neck rupture, that is, at 
the exit point of the fission valley, the mass distribu- 
tion in the spontaneous fission of ^^^Fm is inferred to 
be asymmetric and peaked at A^ w 116 and A^ ~ 140, 
only one mass unit away from the experimental value for 
the heavy fragment Ah « 141 [J. This property has 
also been successfully described within the macroscopic- 
microscopic FRLDM model by MoUer et al. 13] who 
found Ah « 140. Moreover, the calculations by Warda 
et al. 16] have shown the same behavior in the transi- 
tion from the ground-state symmetric path to the asym- 
metric fission valley. This can be seen in their Fig. 5 
where their cuts in the Qn direction at various elonga- 
tions Q2 = Q20/2 are similar to the HFBCS cuts along 
Qso, since Qn and Qao turn out to be in a one-to-one 
correspondence along the cuts of Fig. Inhere. It is worth 
mentioning that the portion of the symmetric path be- 



tween Q2 ~ 90 and Q2 ~ 130 b of their Fig. 3 does not 
correspond to local minima in the Qn direction, as can 
be seen in their Fig. 5. Although the authors of Ref. [Ig 
found a clear left-right reflection asymmetry in their nu- 
clear shapes along the EF path, they reported that the 
two nascent fragments have nearly equal masses, which 
they inferred in an unclear way from the integrated par- 
ticle number as a function of z plotted in their Fig. 4b. 
On the contrary, assuming that the neck most likely rup- 
tures where its radius is minimum, the Figs. 4a and 4b 
of Ref. [la] rather indicate that the heavy fragment has 
a mass of about 136 and an atomic number of about 52, 
in a much better agreement with the HFBCS results and 
the experimental data jj]. 

Contrary to ^^^Fm, the most favorable exit channel 
(i.e, the lowest continuous path) beyond the superde- 
formed minimum in ^^*Fm goes along the only sym- 
metric path which eventually forks, instead of following 
the asymmetric EF valley because of the ridge separat- 
ing them (at least 1.5 MeV high according to Warda 
et al. Id]). Whichever valley is eventually followed by 
the spontaneously fissioning nucleus towards scission, the 
outcome is the same in terms of mass fragmentation since 
in either valley the configurations are left-right symmet- 
ric. This leads thus to a symmetric mass distribution, 
as was obtained in all the other theoretical sudies and in 
agreement with experiment ja, 13 • 



10 



B. Bimodal fission in ^^**Fm 

Let us now relate the properties of the valleys obtained 
in the potential-energy surface of ^^^Fm to the fragment 
total-kinetic-energy and mass distributions of this iso- 
tope. 

The two symmetric valleys present a major difference 
associated with the nuclear shapes. Whereas the fission- 
ing nucleus develops already in the early stages of the 
CF valley a narrow neck connecting two nearly spherical 
nascent fragments (see Fig. I10|l . the neck of a fission- 
ing nucleus following the symmetric EF valley persists 
over a much wider range of total elongation Q2Q, with 
very elongated nascent fragments (see Fig. Ill|) . From a 
geometrical argument and approximating the total frag- 
ment kinetic energy by the Coulomb interaction energy 
(the dominant contribution) at scission, we can deduce 
that the fragments formed in the descent of the fission- 
ing nucleus along the symmetric CF valley have a much 
higher kinetic energy than those associated with fission 
events from the symmetric EF valley. In both cases the 
fragment-mass distribution is expected to be symmetric 
since the parent nucleus fissions into two identical frag- 
ments in both valleys. More specifically, in the case where 
258pj^ undergoes fission through compact shapes like the 
ones displayed in Fig. IIUI the high stiffness of the CF 
valley in the mass asymmetry direction, resulting from 
the strong shell effects in the nearly spherical and dou- 
bly magic nascent fragments, produces a mass distribu- 
tion much narrower than the one corresponding to fission 
events through the symmetric EF valley. 

These arguments show that the two HFBCS symmetric 
valleys are consistent with the kinetic-energy and mass- 
distribution properties of the two modes experimentally 
observed. Therefore it is natural to identify the symmet- 
ric EF path with the low-TKE mode, and the symmetric 
CF path with the high-TKE mode. Another argument 
in favor of the interpretation of the symmetric EF val- 
ley relies on the similarity between the nuclear shapes 
of Fig. ^] and those obtained in the liquid-drop model, 
with which the experimentally measured properties of the 
low-energy mode are consistent . 

To provide more quantitative grounds to this interpre- 
tation, a discussion is now devoted to the estimation of 
the fragment total kinetic energy associated with each 
fission path. Let us assume that TKE is given by 



-Dc.m. = 13.0 fm, Q20 = 226.02 b 



TKE = E, 



(int) 
C 



e: 



(int) 



E- 



(sc) 



where a small nuclear interaction energy E 



(int) 
p(sc) 



(25) 
and 



a somewhat larger prescission kinetic energy Ej^ ' are 
added to the dominant Coulomb interaction energy 



E. 



(int) 



These three contributions have to be calculated 



just after neck rupture, that is, at scission. The Coulomb 
part does not pose any particular problem, so I focus on 
how to calculate the prescission kinetic and nuclear con- 
tributions. The relative energy of a scission configuration 
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FIG. 10: Sequence of nuclear shapes along the symmetric CF 
valley of ^^Spm between Dc.m. = 13.0 fm and Dc.m. = 15.0 fm. 
The solid lines correspond to equal nuclear density contours 
ranging from p(r) — 0.01 fm~^ (outermost contour) to p(r) = 
0.15 fm~ (innermost contour) with 0.01 fm~ steps. The 
thick solid line marks the density contour at half the satu- 
ration density, namely p(r) — psat/2 = 0.08 fm"''. 



Ai^sc with respect to the initial energy of the fissioning 
nucleus (its ground-state energy in the present case of 
spontaneous fission) represents the available energy at 
scission. This energy can be shared among collective de- 
grees of freedom other than deformation, essentially ki- 
netic energy in the fission direction E^ , and internal de- 



grees of freedom as an "internal" excitation energy E\. 



(sc) 
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FIG. 11: Same as Fig. llOl for nuclear shapes along the symmet- 
ric EF valley between Dc.m. = 12.5 fm and Dc.m. = 22.0 fm. 



We can thus write: 



-AK, 



_ p(sc) 



E^^> + E, 



(sc) 



(26) 



assuming that a physical scission configuration lies below 
the ground state to be accessible by tunneling from the 
ground-state well {AEgc < 0). Since a static study is 
unable to provide the actual partitioning, I postulate an 
equipartition of —AEsc/'2 between its two components 



E, 



(sc) 



E, 



(sc) 



A^.„ 



(27) 



This approximation leads to values of prescission kinetic 
energy of about 10 MeV for the scission configurations 



considered in Tab. lIIII This is of the same order as the one 
assumed by Brosa et al. [35| and consistent with the ex- 
perimental average TKE-values, as well as the one found 
by Bonasera |36j| in semi-classical dynamical calculations 
and by Abe et al. |37| | from dynamical calculations based 
on the two-dimensional Langevin equation and the one- 
body dissipation mechanism. As for the nuclear inter- 
action between the two fragments, it can be in principle 
calculated with the Skyrme force. This requires one to 
disentangle the three contributions of the nuclear energy 
of the whole system, namely the self-energies of the two 
fragments and the interaction energy between the frag- 
ments. This is ambiguous since one has to unfold the 
local densities into two sets of densities localized each in 
one of the fragments (with some overlap in the neck). 
However Pomorski and Dietrich [Sg did the calculation 
in the case of two spherical nuclei and showed that the 
resulting potential is similar to the folded Yukawa-plus- 
exponenti al ( YPE) potential proposed by Krappe, Nix 
and Sierk J33 • For this reason I use the latter potential, 
with the parameters of Ref . |43 ■ 

The actual calculation of the different contributions of 
TKE for the three paths requires to determine a scission 
configuration for each path. As for the symmetric CF 
valley, I postulate that the shape at Dc.m. = 15 fm in 
Fig-Elis the scission-point configuration. This is consis- 
tent with the criterion of Goutte et al. |2J| that, at scis- 
sion, the nuclear density in the neck, a,t z = Zneck defined 
in Sect. 2, is 0.01 fm^'^. In the CF valley this corresponds 
to Dc.m. ~ 15.0 fm. In fact it is remarkable that a scis- 
sion point lies in the bottom of a valley. In contrast, no 
scission configurations are found along the symmetric or 
the asymmetric EF paths. It becomes more difficult and 
ambiguous to assign a scission point to each of these EF 
paths. However, in order to obtain an estimate of the 
associated kinetic energies, a kind of "sudden approxi- 
mation" is used for an approximate calculation of i?Q . 
The two nascent fragments are approximated at the exit 
point by the equivalent coaxial spheroids having the same 
elongations and root-mean-square radii as the actual frag- 
ments. They result from a sudden neck rupture assumed 
to preserve the mass and charge asymmetries as well as 
the center-of-mass distance. The Coulomb interaction 
energy is calculated from the exact analytical expression 
of Quentin 41], whereas the nuclear interaction energy is 
calculated with the YPE potential by numerical integra- 
tion. 



The calculated results of E'^'"* , i?^"* i and E^'^' for 
each path are given in Tab. IIIII together with the re- 
sulting kinetic-energy values rounded to the nearest in- 
teger. The TKE-values for the symmetric CF and EF 
valleys lie in the ranges of the experimental high- and 
low-energy modes, respectively (see Fig. 7 of Ref. |2j). 
As for the asymmetric EF path, the HFBCS total ki- 
netic energy is less than 220 MeV, which is consistent 
with the conclusion from Fig. 8 of Ref. |3| that the frag- 
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TABLE III; Characteristics of the approximate scission con- 
figurations for each fission path: center-of-mass distance Dc.m. 
in fm, corresponding total elongation Q20 in barns, deforma- 
tion parameter (3 (dimensionless) defined by Zhao et al. |4^ . 
Coulomb interaction energy Eq , nuclear interaction energy 
E^'^ and prescission kinetic energy E^'^ in MeV, together 
with their sum TKE and the most-probable heavy-fragment 
mass Ah- 
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Dc.u,. 


Q20 


13 


p(int) 


„(int) 

E-f^ 


E^-^ 


TKE Ah 


sym. OF 


15.0 


289.26 


1.27 


238.6 


-1.7 


7.1 


244 129 


asym. EF 


18.2 


462.14 


1.55 


204.7 


-1.4 


11.7 


215 141 


sym. EF 


22.0 


714.36 


1.86 


178.5 


-5.2 


10.4 


184 129 



ments associated with the fission events above Ah = 140 
in the mass distribution have a TKE-value lower than 
220 MeV. The authors of Ref. |9|| also reported that 
the mass-yield curve obtained by selecting spontaneous- 
fission events with TKE < 200 MeV becomes asymmet- 
ric, from which it can be inferred that a significant num- 
ber of symmetric pairs of fragments have a TKE-value 
between 200 and 220 MeV. However there is no evidence 
showing which type of symmetric configurations (elon- 
gated or compact) dominates in this kinetic energy range, 
or what the average TKE-value of fragment pairs with 
Ah ^ 140 is. Therefore, based on kinetic-energy consid- 
erations, it seems possible that the asymmetric EF valley 
contributes to feed the low-energy mode, but probably 
less so than the symmetric EF one because of the ridge 
separating them. 

Two recent studies support this suggestion. On the 
one hand, Zhao et al. |42. |43 | deduced a deformation 
parameter of the scission configurations /?, which gives 
a measure of the deviation from two touching spheres, 
from experimental average total-kinetic-energy system- 
atics. The value f3 = 1.49 that these authors obtained 
from the average TKE-value of the low-energy mode in 
258pj^ is close to the average value /3asym — 1-53 ± 0.02 
corresponding to the asymmetric mode throughout the 
actinide region. It is interesting to note that the value 
I obtained for the asymmetric EF path /3hfbcs — 1-55 
is compatible with this systematics. On the other hand, 
Asano et al. [20j performed more recently dynamical cal- 
culations of the fragment kinetic-energy and mass dis- 
tributions of ^^"'Fm as well as the fragment-mass dis- 
tributions of ^'^^Fm and ^^*Fm, at an excitation en- 
ergy of the compound nucleus of 10 MeV. These au- 
thors found that their distributions can be decomposed 
into three modes: a mass symmetric high-energy mode 
(TKE = 232.1 MeV), a mass asymmetric, low-energy 
mode (A^ = 147.0, TKE = 200.8 MeV) and a symmet- 
ric, very low-energy mode (TKE = 171.7 MeV). Their 
calculated deformation parameters of the scission con- 



figurations associated with the first two modes are in 
rood agreement with the systematics of Zhao et al. 
I43 |. Despite the results obtained by Asano and col- 
laborators correspond to a different isotope at a higher 
compound-nucleus excitation energy, they can be consid- 
ered similar to those for ^^^Fm reported in Tab. IIIII 

Contrary to the above interpretation proposed for the 
calculated fission paths, Warda and collaborators |la.ll7| 
identified the low-kinetic-energy mode with the asymmet- 
ric EF path only, since they do not seem to have found 
a symmetric EF path. In the same way as for ^^^Fm, 
they accounted for the symmetric character of the corre- 
sponding mass distribution by a mass symmetric division 
associated with left-right reflection asymmetric shapes. 
This explanation is not supported by the present HF- 
BCS calculations. However, it would be very interesting 
to compare the exit points of the asymmetric EF valleys 
obtained in both models. Indeed the finite-range effect 
of the Gogny effective force may play a role in the nu- 
clear interaction energy, therefore impacting the scission 
configurations and the total kinetic energies. 

Finally the similar abundance experimentally observed 
for the two fission modes still needs to be explained. Since 
a static model cannot predict the branching ratio, I can 
only provide the following plausible qualitative argument. 
As we can see in Fig. [7| the symmetric EF path appears 
between Dc.m. = 12.0 fm and Dc.m. = 12.5 fm, at about 
the same energy as the CF path. This seems to indicate 
an equally important feeding of both valleys, leading to 
an expected branching ratio of about 1. 



V. CONCLUSION 

The features of the potential-energy surface of 
the ^^^Fm and ^^^Fm isotopes calculated within the 
HF(SkM*)-fBCS(G) model successfully account for 
the experimentally observed asymmetric-to-symmetric 
transition in the mass distribution for spontaneous 
fission as well as most of the measured properties of 
the bimodal spontaneous fission in ^^^Fm. The HFBCS 
results suggest a different interpretation from the one 
proposed by Warda et al. 16j for the low-energy mode in 
the spontaneous fission of ^^^Fm. This mode seems to be 
better understood as a combination of a dominant one 
corresponding to symmetric very elongated scission con- 
figurations, and another one associated with asymmetric 
rather elongated scission shapes (with a nearly spherical 
heavy fragment). The estimates of the total kinetic 
energy for each component are compatible with the 
experimental data |9| as well as with recent dynamical 
calculations [2y|. As an alternate dynamical approach, 
the virial-theorem-based approach, developped and ini- 
tially applied to heavy-ion collisions by I. N. Mikhailov 
and collaborators ^^, can also be applied to fission. 
This work is underway for the symmetric fission of the 
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258pj^ isotope with microscopically calculated ingredi- 
ents, namely the potential-energy surface obtained in 
the present study and the inertia parameters calculated 
in the HFBCS framework as in Rcf. lisll. 
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